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We describe a novel oscillator that converts continuous light energy into sta- 
ble and spectrally pure microwave signals. This light-induced microwave oscillator 
(LIMO) consists of a pump laser and a feedback circuit, including an intensity 
modulator, an optical fiber delay line, a photodetector, an amplifier, and a filter . 
We develop a quasilinear theory and obtain expressions for the threshold condition, 
the amplitude, the frequency ; the line width, and the spectral power density of 
the oscillation. We also present experimental data to compare with the theoret- 
ical results. Our findings indicate that the LIMO can generate ultrastable, spec- 
trally pure microwave reference signals up to 75 GHz with a phase noise lower than 
— 140 dBc/Hz at 10 kHz. 


I. Introduction 

Oscillators are devices that convert energy from a continuous source to a periodically varying signal. 
They represent the physical realization of a fundamental basis of all physics, the harmonic oscillator, and 
they are perhaps the most widely used devices in modern day society Today a variety of oscillators — 
mechanical [1] (such as the pendulum), electromagnetic (such as LC circuit [2,3] and cavity-based [4]), and 
atomic (such as masers [5] and lasers [6]) — provides a diverse range in the approximation to the realization 
of the ideal harmonic oscillator. The degree of the spectral purity and stability of the output signal of the 
oscillator is the measure of the accuracy of this approximation and is fundamentally dependent on the 
energy storage ability of the oscillator, determined by the resistive loss (generally frequency dependent) 
of the various elements in the oscillator. 


An important type of oscillator widely used today is the electronic oscillator. The first such oscillator 
was invented by L. De Forest [2] in 1912, shortly after the development of the vacuum tube. In this 
triode-based device known as the van der Pol oscillator [3], the flux of electrons emitted by the cathode 
and flowing to the anode is modulated by the potential on the intervening grid. This potential is derived 
from the feedback of the current in the anode circuit containing an energy storage element (i.e., the 
frequency-selecting LC filter) to the grid, as shown in Fig. 1(a). Today the solid-state counterparts 
of these valve oscillators based on transistors are pervasive in virtually every application of electronic 
devices, instruments, and systems. Despite their widespread use, electronic oscillators, whether of the 
vacuum-tube or the solid-state variety, are relatively noisy and lack adequate stability for applications 
where very high stability and spectral purity are required. The limitation to the performance of electronic 
oscillators is due to ohmic and dispersive losses in various elements in the oscillator, including the LC 
resonant circuit. 
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Fig. 1. Comparison of two types of oscillators: (a) van der Pol oscillator and (b) light-induced microwave oscillator. 


For approximately the past 50 years, the practice of reducing the noise in the electronic oscillator by 
combining it with a high-quality factor (Q) resonator has been followed to achieve improved stability and 
spectral purity. The Q is a figure of merit for the resonator given by Q — 2n fr^ where is the energy 
decay time that measures the energy storage ability of the resonator and / is the resonant frequency. 
High-Q resonators used for stabilization of the electronic oscillator include mechanical resonators, such 
as quartz crystals [7,8]; electromagnetic resonators, such as dielectric cavities [9]; and acoustic [10] and 
electrical delay lines, where the delay time is equivalent to the energy decay time, r^, and determines 
the achievable Q. This combination with a resonator results in hybrid-type oscillators referred to as 
electromechanical, electromagnetic, or electro-acoustic, depending on the particular resonator used with 
the oscillator circuit. The choice of the particular resonator is generally determined by a variety of factors, 
but for the highest achievable Q’s at room temperatures, the crystal quartz is the resonator of choice for 
the stabilization of the electronic oscillator. However, because quartz resonators have only a few high-Q 
resonant modes at low frequencies [7,8], they have a limited range of frequency tunability and cannot be 
used to directly generate high-frequency signals. 


In this article, we introduce a novel photonic oscillator [11,12] characterized by spectral purity and 
frequency stability rivaling the best crystal oscillators. This oscillator, shown schematically in Fig. 1(b), 
is based on converting the continuous light energy from a pump laser to radio frequency (RF) and 
microwave signals, and thus we refer to it with the acronym LIMO, for the light-induced microwave 
oscillator. The LIMO is fundamentally similar to the van der Pol oscillator, with photons replacing 
the function of electrons, an electro-optic (E/O) modulator replacing the function of the grid, and a 
photodetector replacing the function of the anode. The energy storage function of the LC circuit in the 
van der Pol oscillator is replaced with a long fiber-optic delay line in the LIMO. 

Despite this close similarity, the LIMO is characterized by significantly lower noise and very high 
stability, as well as by other functional characteristics that are not achieved with the electronic oscillator. 
The superior performance of the LIMO results from the use of electro-optic and photonic components that 
are generally characterized by high efficiency, high speed, and low dispersion in the microwave frequency 
regime. Specifically, currently there are photodetectors available with quantum efficiency as high as 
DO percent that can respond to signals with frequencies as high as 110 GHz [13]. Similarly, E/O modulators 
with a 75-GHz frequency response are also available [14]. Finally, the commercially available optical fiber, 
which has a small loss of 0.2 dB/km for 1550 nm light, allows long storage time of the optical energy 
with negligible dispersive loss (loss dependent on frequency) for the intensity modulations at microwave 
frequencies. 
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The LIMO may also be considered as a hybrid oscillator in so far as its operation involves both light 
energy and microwave signals. Nevertheless, as a hybrid oscillator, the LIMO is unique in that its output 
may be obtained both directly as a microwave signal or as intensity modulation of a optical carrier. This 
property of the LIMO is quite important for applications involving optical elements, devices, or systems 

in]. 

The ring configuration, consisting of an electro-optic modulator that is fed back with a signal from the 
detected light at its output, has been previously studied by a number of investigators interested in the 
nonlinear dynamics of bistable optical devices [15-19]. The use of this configuration as a possible oscillator 
was first suggested by A. Neyer and E. Voges [20]. Their investigations, however, were primarily focused on 
the nonlinear regime and the chaotic dynamics of the oscillator. This same interest persisted in the work of 
T. Aida and P. Davis [21], who used a fiber wave guide as a delay line in the loop. Our studies, by contrast, 
are specifically focused on the stable oscillation dynamics and the noise properties of the oscillator. The 
sustainable quasilinear dynamics, both in our theoretical and experimental demonstrations, are arrived 
at by the inclusion of a filter in the feedback loop to eliminate harmonics generated by the nonlinear 
response of the E/O modulator. This approach yields stable, low-noise oscillations that closely support 
the analytical formulation presented here. 

In this article, we first describe the oscillator and identify the physical basis for its operation. We then 
develop a quasilinear theory for the oscillator dynamics and the oscillator noise. Results of the theory 
are then compared with experimental results. 


il. Description of the Oscillator 

The LIMO utilizes the transmission characteristics of a modulator together with a fiber-optic delay 
line to convert light energy into stable, spectrally pure RF /microwave reference signals. A detailed view 
of the construction of the oscillator is shown schematically in Fig. 2. In this depiction, light from a laser 
is introduced into an E/O modulator, the output of which is passed through a long optical fiber, and 
detected with a photodetector. The output of the photodetector is amplified and filtered and fed back to 
the electric port of the modulator. This configuration supports self-sustained oscillations at a frequency 
determined by the fiber delay length, bias setting of the modulator, and the bandpass characteristics of 
the filter. It also provides for both electrical and optical outputs, a feature which would be of considerable 
advantage to photonics applications. 

We use a regenerative feedback model to analyze the spectral properties of the LIMO. Similar methods 
have been successfully used to analyze lasers [5] and surface acoustic wave oscillators [22]. The conditions 
for self sustained oscillations include coherent addition of partial waves each way around the loop and a 
loop gain exceeding losses for the circulating waves in the loop. The first condition implies that all signals 
that differ in phase by some multiple of 2n from the fundamental signal may be sustained. Thus, the 
oscillation frequency is limited only by the characteristic frequency response of the modulator and the 
setting of the filter, which eliminates all other sustainable oscillations. The second condition implies that 
with adequate light input power, self sustained oscillations may be obtained, without the need for the 
RF/microwave amplifier in the loop. These characteristics, which are expected based on the qualitative 
analysis of the oscillator dynamics, are mathematically derived in the following sections. 


III. Quasilinear Theory of the LIMO 

In the following sections, we introduce a quasilinear theory to study the dynamics and noise of a 
LIMO. In the discussion, we assume that the E/O modulator in the oscillator is of the Mach-Zehnder 
type. However, the analysis of oscillators made with different E/O modulators may follow the same 
procedure. The flow of the theory is as follows: First, the open-loop characteristics of a photonic link 
consisting of a laser, a modulator, a fiber delay, and a photodetector are determined. We then close the 
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Fig. 2. Detailed construction of a LIMO. Optical injection and RF injection ports are 
supplied for synchronizing the oscillator with an external reference by either optical 
injection locking or electrical injection locking [12]. The bias port and the fiber stretcher 
can be used to fine tune the oscillation frequency [12]. Noise in the oscillator can be 
viewed as being injected from the input of the amplifier. 


loop back into the modulator and invoke a quasilinear analysis by including a filter in the loop. This 
approach leads us to a formulation for the amplitude and the frequency of the oscillation. In the next 
step, we consider the influence of the noise in the oscillator, again assisted by the presence of the filter, 
which limits the number of circulating Fourier components. We finally arrive at an expression for the 
spectral density of the LIMO that would be suitable for experimental investigations. 

A. Oscillation Threshold 

The optical power from the E/O modulator's output port that forms the loop is related to an applied 
voltage V in (t) by 


P(t) = 



Tf sill 7 T 


V<n(0 

K 



( 1 ) 


where a is the fractional insertion loss of the modulator, V n is its half-wave voltage, Vb is its bias voltage, 
P 0 is the input optical power, and 77 determines the extinction ratio of the modulator by (1 4- 77) /(I - 77). 

If the optical signal P(t) is converted to an electric signal by a photodetector, the output electric signal 
after an RF amplifier is 


Vout(t) = pP(t)RG A = V p h{ 1 - 77 sin 7r 


Vin(t) Vb 

Vk V ; 


( 2 ) 


where p is the responsivity of the detector, R is the load impedance of the photodetector, G A is the 
amplifier’s voltage gain, and V ph is the photovoltage, defined as 


v ph = RG a = I vh RG A 


( 3 ) 


with I ph = aP 0 p /2 as the photocurrent. The LIMO is formed by feeding the signal of Eq. (2) back to the 
RF input port of the E/O modulator. Therefore, the small signal open- loop gain, Gs, of the LIMO is 
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V; 7l =0 


( 4 ) 


G s = 


dVput 

~dV~ 



The highest small signal gain is obtained when the modulator is biased at quadrature, that is, when 
Vb = 0 or Vn- From Eq. (4), one may see that Gs can be either positive or negative, depending on 
the bias voltage. The modulator is said to be positively biased if Gs > 0; otherwise it is negatively 
biased. Therefore, when V B = 0, the modulator is biased at negative quadrature, while when V B = V n , 
the modulator is biased at positive quadrature. Note that in most externally modulated photonic links 
the E/O modulators can be biased at either positive or negative quadrature without affecting their 
performance. However, as will be seen next, the biasing polarity will have an important effect on the 
operation of the LIMO. 

In order for the LIMO to oscillate, the magnitude of the small signal open-loop gain must be larger 
than unity. From Eq. (4), we immediately obtain the oscillation threshold of the LIMO to be 

Kr 

Vph 7TT7|cos(?rV B /%)| ^ 

For the ideal case in which 77 = 1 and V B — 0 or V n , Eq. (5) becomes 

Vph = — (6) 

7 r 

It is important to notice from Eqs. (3) and (6) that the amplifier in the loop is not a necessary condition 
for oscillation. So long as I p hR > / 7r is satisfied, no amplifier is needed ( Ga = 1)- It is the optical 

power from the pump laser that actually supplies the necessary energy for the photonic oscillator. This 
property is of practical significance because it enables the LIMO to be powered remotely using an optical 
fiber. Perhaps more significantly, however, the elimination of the amplifier in the loop also eliminates 
the amplifier noise, resulting in a more stable oscillator. For a modulator with a V n of 3.14 V and 
an impedance, R , of 50 a photocurrent of 20 mA is required for sustaining the photonic oscillation 
without an amplifier. This corresponds to an optical power of 25 mW, assuming the responsivity, p, of 
the photodetector to be 0.8 A/W. 

B. Linearization Of the E/O Modulator’s Response Function 

In general, Eq. (2) is nonlinear. If the electrical input signal Vi n {t) to the modulator is a sinusoidal 
wave with an angular frequency of u, an amplitude of V 0 , and an initial phase of /?, 


V in (t) = V 0 sm(ust-h/3) (7) 

then the output at the photodetector, V out (t ) , can be obtained by substituting Eq. (7) in Eq. (2) and 
expanding the left-hand side of Eq. (2) with Bessel functions: 


Fout(^) — Vph 



cos(2 mut 4- 2m/?) 


- 2r] cos 




sin [(2m + 1 )ujt + (2m + 1)/?] 


( 8 ) 
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It is clear from Eq. (8) that the output contains many harmonic components of lj. 

The output can be linearized if it passes through an RF filter with a bandwidth sufficiently narrow to 
block all harmonic components. The linearized output can be obtained easily from Eq. (8): 


Vout(t) = G(V 0 )V in (t) 


where the voltage gain coefficient, G(Vq), is defined as 


2V 

G(V 0 ) = G s —^J l 

7T Vq 



(9) 


(10a) 


It can be seen that the voltage gain, G(V 0 ), is a nonlinear function of the input amplitude, V Q , and its 
magnitude decreases monotonically with V Q . However, for a small enough input signal (V Q << V n and 
J^TiVo/V-n) - 7vV 0 /2V n ), we can recover from Eq. (10a) the small signal gain: G(V 0 ) — Gs- 

If we expand the left-hand side of Eq. (2) with Taylor series, the gain coefficient can be obtained as 


G(V 0 ) = G s 



(ttVo 

\2Vn 




(10b) 


It should be kept in mind that, in general, G ( V () ) is also a function of the frequency, u, of the input 
signal, because V p h is linearly proportional to the gain of the RF amplifier and the responsivity of the 
photodetector, which are all frequency dependent. In addition, the V n of the modulator is also a function 
of the input RF frequency. Furthermore, the frequency response of the RF filter in the loop can also be 
lumped into G{V 0 ). In the discussions below, we will introduce a unitless complex filter function, F(cj), 
to explicitly account for the combined effect of all frequency-dependent components in the loop while 
treating G(V 0 ) to be frequency independent: 

F{u) = F(u>)e t<t>{u] (11) 

where 0(u;) is the frequency- dependent phase caused by the dispersive component in the loop and F(&) 
is the real normalized transmission function. Now Eq. (9) can be rewritten in complex form as 


V out (t) ^ F(u;)G(V 0 )V in (^t) 


(12) 


where V in and V out are complex input and output voltages. Note that although Eq. (12) is linear, the 
nonlinear effect of the modulator is not lost — it is contained in the nonlinear gain coefficient, G(V a ). 


C. Oscillation Frequency and Amplitude 

In this section, we derive the expressions for the amplitude and frequency of the LIMO. Like other 
oscillators, the oscillation of a LIMO starts from noise transient, which is then built up and sustained 
with feedback at the level of the oscillator output signal. We derive the amplitude of the oscillating 
signal by considering this process mathematically. The noise transient can be viewed as a collection of 
sine waves with random phases and amplitudes. To simplify our derivation, we use this noise input with 
the linearized expression of Eq. (12) for the loop response. Because Eq. (12) is linear, the superposition 
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principle holds, and we can analyze the response of the LIMO by first inspecting the influence of a 
single-frequency component of the noise spectrum: 


V in (w,t) = V in (ui)e iuJt 


(13) 


where V in (u)) is a complex amplitude of the frequency component. 

Once the noise component of Eq. (13) is in the oscillator, it would circulate in the loop, and the 
recurrence relation of the fields from Eq. (12) is 


Vnfat) = F(u)G{V 0 )V n ^t^ r') 


(14) 


where r' is the time delay resulting from the physical length of the feedback and n is the number of times 
the field has circulated around the loop, with V n= o(cj,t) = V in (u>,t). In Eq. (14), the argument V Q in 
G(V 0 ) is the amplitude of the total field (the sum of all circulating fields) in the loop. 

The total field at any instant of time is the summation of all circulating fields. Therefore, with the 
input of Eq. (13) injected in the oscillator, the signal measured at the RF input to the modulator for the 
case when the open-loop gain is less than unity can be expressed as 


V(w,t) = G a V in (u>) Y, [F(u)G(V 0 )] n e iu{t - nT,) 

71 = 0 


GgVin^Y 

1 - F(u)G(V 0 )e-^ T ' 


(15) 


For a loop gain below threshold and with a small V a , G(V 0 ) is essentially the small signal gain, Gs , given 
by Eq. (4). 

The corresponding RF power of the circulating noise at frequency to is, therefore, 

|VW)| a <%|VinMI 2 /(2fl) , 1(n 

(u} 2 R 1 + \F(w)G(V 0 )\ 2 - 2F’(w)|G(C,)| cos[u;r' + <t>(u>) + 4>o}' K 1 


where 0 O = 0 if G(V 0 ) > 0 and (j) 0 = r if G{V 0 ) < 0. 

For a constant V in ( cj), the frequency response of a LIMO has equally spaced peaks similar to those of 
a Fabry-Perot resonator, as shown in Fig. 3. These peaks are located at the frequencies determined by 


u^t' + (p{^k) + 4>o ~ 2/c7t A; = 0, 1,2,--- (17) 

where k is the mode number. In Fig. 3, each peak corresponds to a frequency component resulting from 
the coherent summation of all circulating fields in the loop at that frequency. As the open-loop gain 
increases, the magnitude of each peak becomes larger and its shape becomes sharper. These peaks are 
the possible oscillation modes of the LIMO. When the open- loop gain is larger than unity, each time a 
noise component at a peak frequency travels around the loop, it is amplified and its amplitude increases 
geometrically — an oscillation is started from noise. 

Because an RF filter is placed in the loop, the gain of only one mode is allowed to be larger than unity, 
thus selecting the mode that is allowed to oscillate. Because of the nonlinearity of the E/O modulator 
or the RF amplifier, the amplitude of the oscillation mode cannot increase indefinitely. As the amplitude 
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Fig. 3. Illustration of the oscillator's output spectra below and above the threshold. 


increases, higher harmonics of the oscillation will be generated by the nonlinear effect of the modulator or 
the amplifier, at the expense of the oscillation power, and these higher harmonics will be filtered out by 
the RF filter. Effectively, the gain of the oscillation mode is decreased according to Eq. (10) until the gain 
is, for all practical measures, equal to unity, and the oscillation is stable. As will be shown later, because 
of the continuous presence of noise, the closed-loop gain of an oscillating mode is actually less than unity 
by a tiny amount on the order of 10 — 10 , which ensures that the summation in Eq. (15) converges. 

In the discussion that follows, only one mode k is allowed to oscillate, and so we will denote the 
oscillation frequency of this mode as f osc or uj 0 sc{u 0 sc — 27 r/ 05C ), its oscillation amplitude as V osc1 and its 
oscillation power as P 0 sc(Posc = Vosc! 2ft)- In this case, the amplitude, Vo, of the total field in Eq. (16) 
is just the oscillation amplitude, V osc , of the oscillating mode. If we choose the transmission peak of the 
filter to be at the oscillation frequency, cj 05C , and so F(lj osc ) — 1, the oscillation amplitude can be solved 
by setting the gain coefficient, |G(\^, SC )|, in Eq. (16) at unity. From Eq. (10a), this leads to 


J i 



1 7 tV osc 

2\G S \ V n 


(18a) 


In deriving Eq. (18a), we have assumed that the RF amplifier in the loop is linear enough that the 
oscillation power is limited by the nonlinearity of the E/O modulator. The amplitude of the oscillation 
can be obtained by solving Eq. (18a) graphically, and the result is shown in Fig. 4(a). Note that this 
result is the same as that obtained by Neyer and Voges [20] using a more complicated approach. 

If we use Eq, (10b), we can obtain the approximated solutions of the oscillation amplitude: 


Vn, 




i 

W\ 


thirdorder expansion (18b) 


v oat 


2V3K ^ 


1 

7s 



1/2 


fifthorder expansion 


(18c) 


The threshold condition of |Gs| > 1 is clearly indicated in Eqs. (18b) and (18c). Figure 4(a) shows 
the normalized oscillation amplitude as a function of |Gs|, obtained from Eqs. (18a), (18b) and (18c), 
respectively. Comparing the three theoretical curves, one can see that, for |Gs| < 1.5, the third-order 
expansion result is a good approximation. For |G$| < 3, the fifth-order expansion result is a good 
approximation. 
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Fig. 4. The normalized oscillation amplitude of a LIMO as a function 
of small signal gainJG^: (a) theoretical calculation using Eqs. (18a), 
(18b), and (18c) and (b) experimental data and curve fitting to Eqs. 
(18b) and (18c). 


The corresponding oscillation frequency, f osc = = Uk/ 2tt, can be obtained from Eq. (17) as 


_ , (* + 1/2) 
- Jk ~ 

T 

for G(V 0SC ) < 0 

(19a) 

= fk=- 

T 

for G(V r osc ) > 0 

(19b) 


where r is the total group delay of the loop, including the physical length delay, r', of the loop and the 
group delay resulting from dispersive components (such as an amplifier) in the loop, and it is given by 


T = T + 


d(f)(uj) 


duj 


( 20 ) 


For all practical purposes, J\{^V osc IV^) > 0 or V osc /Vn < 1.21, and the sign of G(V osc ) is determined by 
the small signal gain, G s . It is interesting to notice from Eqs. (19a) and (19b) that the oscillation frequency 
depends on the biasing polarity of the modulator. For negative biasing (G s < 0), the fundamental 
frequency is l/(2r), while for positive biasing (G s > 0), the fundamental frequency is doubled to 1/r. 
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D. The Spectrum 

The fundamental noise in a LIMO consists of the thermal noise, the shot noise, and the laser’s intensity 
noise, which for the purpose of analysis can be viewed as all originating from the photodetector. Since 
the photodetector is directly connected to the amplifier, the noise can be viewed as entering the oscillator 
at the input of the amplifier, as shown in Fig. 2. 

We compute the spectrum of the oscillator signal by determining the power spectral density of noise 
in the oscillator. Let Pn{u)) be the power density of the input noise at frequency <j; we have 


p N (uj)Af = 


lUnMI 2 

27? 


( 21 ) 


where Af is the frequency bandwidth. Substituting Eq. (21) in Eq. (16) and letting F(lj osc ) = 1, we 
obtain the power spectral density of the oscillating mode k: 


SrfU') 


P(f ) 

A fPosc 


PnG'a/ P osc 

1 + \F(f')G(V osc )\i - 2F(f')\G(V osc )\cos(27Tf'T) 


( 22 ) 


where = — uj osc ) /2tx is the frequency offset from the oscillation peak, / 05C « In deriving Eq. (22), 

both Eqs. (17) and (20) are used. 

By using the normalization condition 


/ oo /*l/2r 

SrfWW « / S RF (f')df = 1 (23) 

-oo «/— l/2r 

we obtain 

1 - |G(K osc )| 2 « 2 [1 - |G(V osc )|] = (24) 

Note that in Eq. (23) we have assumed that the spectral width of the oscillating mode is much smaller 
than the mode spacing, 1/r, of the oscillator, so that the integration over 1/r is sufficiently accurate. In 
addition, in the derivation, we have assumed that |F(/')| ~ 1 in the frequency band of integration. 

Typically, pjy ~ 10“ 17 mW/Hz, P osc ~ 10 mW, G\ ~ 100, and r ~ 10 -G s. From Eq. (24), one can 
see that the closed-loop gain, IG^V^c)!, of the oscillating mode is less than unity by an amount of 10“ 10 . 
Therefore, the equation 1(2(1^^) | — 1 is sufficiently accurate for calculating the oscillation amplitude, 
Vosc , as in Eqs. (18a), (18b), and (18c). 

Finally, substituting Ecp (24) in Eq. (22), we obtain the RF spectral density of the LIMO: 


SiifU') 


6 

(2 - 6/r) - 2yJ\ - 6/t cos(2nf'T) 


(25) 


where 6 is defined as 
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( 26 ) 


3 _ PnG\ 

Pose 

As mentioned before, p^ is the equivalent input noise density injected into the oscillator from the input 
port of the amplifier, and P osc jG 2 A is the total oscillating power measured before the amplifier. Therefore, 
6 is the input noise- to-signal ratio to the oscillator. 

For the case where 2 7t/'t << 1, we can simplify Eq. (25) by expanding the cosine function in Taylor 
series: 


SrfU') 


6 

(< 5 / 2 t ) 2 + ( 2 t ) 2 ( t /') 2 


(27) 


Equation (27) is a good approximation even for 27 r/'r = 0.7, at which value the error resulting from 
neglecting the higher-order terms in Taylor expansion is less than 1 percent. It can be seen from Eq. (27) 
that the spectral density of the oscillating mode is a Lorentzian function of frequency. Its full width at 
half maximum (FWHM), A /fwhm, is 


A Ifwhm = 


27 r r 2 


_L G aPn 

27 r t 2 P osc 


(28a) 


It is evident from Eq. (28a) that A Jfwhm is inversely proportional to the square of the loop delay time 
and linearly proportional to the input noise- to- signal ratio, 6. For a typical <5 of 10 -16 /Hz and a loop 
delay of 100 ns (20 m), the resulting spectral width is submilli-Hertz. The fractional power contained in 
A Ifwhm is A /fwhmS'/jf(0) = $4 percent. 

From Eq. (28a), one can see that for fixed px and Ga, the spectral width of a LIMO is inversely pro- 
portional to the oscillation power, similar to the famous Schawlow-Townes formula [23,35] for describing 
the spectral width, A viaser, of a laser: 


A l^laser — 


_ Ps 
27T ^~laser^ aser 


(28b) 


where p s = kv is the spontaneous emission noise density of the laser, Pi aS er is the laser oscillation power, 
and Ti aser is the decay time of the laser cavity. However, because, as will be shown in Section III.E, 
both P osc and px are functions of the photocurrent, the statement that the spectral width of a LIMO is 
inversely proportional to the oscillation power is valid only when thermal noise dominates in the oscillator 
at low photocurrent levels. 

From Eq. (28a), the quality factor, Q, of the oscillator is 


Q = 


fo. 


A Ifwhm 



where Qo is the quality factor of the loop delay line and is defined as 


(29) 


Qd = 2 t rf osc T 


(30) 


From Eq. (27), we easily obtain 
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(31a) 


SrfW) 


4 r 2 

6 


i/'l « 


AIfwhm 

2 


SflF(/0 


6 

( 27 t ) 2 ( t /') 2 


l/'l » 


AJfwhm 

2 


(31b) 


It can be shown [24] that for an oscillator with a phase fluctuation much less than unity, its power spectral 
density is equal to the sum of the single side-band phase noise density and the single side-band amplitude 
noise density. In most cases in which the amplitude fluctuation is much less than the phase fluctuation, 
the power spectral density is just the single side-band phase noise. Therefore, it is evident from Eq. (31b) 
that the phase noise of the LIMO decreases quadratically with the frequency offset, /'. For a fixed 
the phase noise decreases quadratically with the loop delay time. The larger the r, the smaller the phase 
noise. However, the phase noise cannot decrease to zero no matter how large r is, because at large enough 
r, Eqs. (27) and (31b) are not valid anymore. From Eq. (27), the minimum phase noise is Sjjj? ~ 6/4 at 
/' = l/2r. For the frequency offset, /', outside of the passband of the loop filter (where F(f f ) = 0), the 
phase noise is simply the noise-to-signal ratio, <5, as can be seen from Eq. (22). 

Equations (25), (27), and (31b) also indicate that the oscillator’s phase noise is independent of the 
oscillation frequency, f osc . This result is significant because it allows the generation of high-frequency and 
low phase- noise signals with the LIMO. The phase noise of a signal generated using frequency-multiplying 
methods generally increases quadratically with the frequency. 

E. The Noise-to-Signal Ratio 

As mentioned before, the total noise density input to the oscillator is the sum of the thermal noise, 
P thermal = 4ksT{NF)\ the shot noise, p s hot — 2 eI p h,R\ and the laser’s relative intensity noise (RIN), 
Prin = N RIN I* h R, densities [25,26]: 


Pn = ikBT(NF) + 2 eI p hR + N R iNl^ h R 


(32) 


where ks is Boltzmann’s constant, T is the ambient temperature, NF is the noise factor of the RE 
amplifier, e is the electron charge, I p h is the photocurrent across the load resistor of the photodetector, 
and Nrin is the RIN of the pump laser. 

From Eqs. (26) and (32), one can see that if the thermal noise is dominant, then 6 is inversely 
proportional to the oscillating power, Pose of the oscillator. In general, P osc is a function of photocurrent, 
Iph, and amplifier gain, G A , as determined by Eqs. (18a), (18b), and (18c), and the noise-to-signal ratio 
from Eq. (26) is thus, 


IG^I 2 4kT(NF) + 2eI ph R 4- N RIN I* h R 
“ 1 - 1/|G S | 4r ? 2 co S 2(7r Vb/V^D^R [ 

111 deriving Eq. (33), Eqs. (4) and (18b) are used. From Eq. (33), one can see that 6 is a nonlinear 
function of the small signal gain of the oscillator. As shown in Fig. 5(a), it reaches the minimum value 
at | G s { — 3/2: 


^min — 


4kT(NF) + 2 eI ph R + N RIN P ph R 
(16/27 )P ph R 


(34) 
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Fig. 5. Calculations of the LIMO's input noise-to-signal ratio (a) as 
a function of small signal gain, IG S I, showing a minimum value at 
lGgl= 1.5 and (b) as a function of photocurrent for different values of 
the laser s RIN noise. 


where rj = 1 and cos^Vjg/V^) = 1 are assumed. The oscillation amplitude at IG 5 I = 3/2 can be obtained 
from Eq. (18b) as 


_ 2\Z2Vtc _ p r 0 T/ 

Kosc — ^ 0.52 Vv 

(?rV3) 


and the corresponding BF power is 


P — 

J O.SC 


4v; 2 _ iop^ dB 

(3tt 2 P) “ 3 


(35a) 


(35b) 


where P^ dB is the input 1 -dB compression power of the E/O modulator [26,27]. From Eq. (35b), one can 
conclude that, in order to have minimum noise, the oscillation power measured at the input of the E/O 
modulator should be 5 dB above the 1 -dB compression power of the modulator. Equation (35a) indicates 
that the noise of the oscillator is at minimum when the oscillating amplitude is roughly half of V n or the 
voltage in the oscillator is varying between the peak and the trough of the sinusoidal transmission curve 
of the E/O modulator. This makes sense because the modulator has its minimum sensitivity to voltage 
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variations at the maximum and the minimum of the transmission curve, and the most likely cause of 
voltage variations in a LIMO is the noise in the loop. 

It is evident from Eq. (34) that the higher the photocurrent, the less the noise-to-signal ratio of the 
oscillator until it flattens out at the laser’s RIN level. Therefore, the ultimate noise-to-signal ratio of 
a LIMO is limited by the pump laser’s RIN. If the RIN of the pump laser in a LIMO is -160 dB/Hz, 
the ultimate noise-to-signal ratio of the oscillator is also —160 dB/Hz, and the signal-to-noise ratio is 
160 dB/Hz. Figure 5(b) shows the noise-to-signal ratio, <5, as a function of photocurrent, / P / M for different 
RIN levels. In the plot, the small signal gain, Gs , is chosen to be a constant of 1.5, which implies that 
when I p f h is increased, the amplifier gain, G A , must be decreased to keep Gs constant. From the figure, 
one can easily see that 8 decreases quadratically with I p h at small I p h and flattens out at the RIN level 
at large I ph . 


F. Effects of Amplifier Nonlinearity 

In the discussions above, we have assumed that the nonlinear distortion of a signal from the E/O 
modulator is more severe than from the amplifier (if any) used in the oscillator, so that the oscillation 
amplitude or power is limited by the nonlinear response of the E/O modulator. Using an engineering 
term, this simply means that the output 1-dB compression power of the amplifier is much larger than the 
input 1-dB compression power of the E/O modulator [26]. 

For cases in which the output 1-dB compression power of the amplifier is less than the input 1-dB 
compression power of the modulator, the nonlinearity of the amplifier will limit the oscillation amplitude, 
V osc (or power P osc ), of the oscillator, resulting in an oscillation amplitude less than that given by 
Eqs. (18a), (18b), and (18c). The exact relation between the oscillation amplitude (or power) and the 
small signal gain, Gs, can be determined using the same linearization procedure as that used for obtaining 
Eqs. (18a), (18b), and (18c) if the nonlinear response function of the amplifier is known. However, all 
the equations in Section III.D for describing the spectrum of the oscillator are still valid, provided that 
the oscillation power in those equations is determined by the nonlinearity of the amplifier. For a high 
enough small signal gain, Gs, the oscillation power is approximately a few dB above the output 1-dB 
compression power of the amplifier. 

In all the experiments below, the output 1-dB compression power of the amplifiers chosen is much 
larger than the input 1-dB compression power of the modulator, so that the oscillation power is limited 
by the modulator. 

IV. Experiments 

A. Amplitude Versus Open-Loop Gain 

We performed measurements to test the level of agreement of the theory described above with exper- 
imental results. In all of our experiments, we used a highly stable diode-pumped Nd:YAG ring laser [35] 
with a built-in RIN reduction circuit [36] to pump the LIMO. The experimental setup for measuring the 
oscillation amplitude as a function of the open-loop gain is shown in Fig. 6(a). Here an RF switch was 
used to open and close the loop. While the loop was open, an RF signal from a signal generator with the 
same frequency as the oscillator was injected into the E/O modulator. The amplitudes of the injected 
signal and the output signal from the loop were measured with an oscilloscope to obtain the open-loop 
gam, which was the ratio of the output amplitude to the injected signal amplitude. The open-loop gain 
was varied by changing the bias voltage of the E/O modulator, by attenuating the optical power of the 
loop, or by using a. variable RF attenuator after the photodetector, as indicated by Eq. (4). When closing 
the loop, the amplitude of the oscillation was conveniently measured using the same oscilloscope. We 
measured the oscillation amplitudes of the LIMO for different open-loop gains at an oscillation frequency 
of 100 MHz, and the data obtained are plotted in Fig. 4(b). It is evident that the experimental data 
agreed well with our theoretical predictions. 
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SIGNAL ANALYZER ANALYZER RF AMPLIFIER 


Fig. 6. Experimental setups for measuring (a) the oscillation amplitude of a LIMO as a function of the 
small signal gain and (b) the phase noise of a LIMO using the frequency discrimination method. 


B. Phase Noise Measurement Setup 

We used a frequency discriminator method [28] to measure the phase noise of the LIMO, and the 
experimental setup is shown in Fig. 6(b). The advantage of this method is that it does not require a 
frequency reference and, hence, can be used to measure an oscillator of any frequency. Using a microwave 
mixer in the experiment, we compared the phase of a signal from the electrical output port of the LIMO 
with its delayed replica from the optical output port. The length of the delay line is important because, 
the longer the delay line, the lower the frequency offset at which the phase noise can be accurately 
measured. On the other hand, if the delay line is too long, the accuracy of the phase noise at a higher 
frequency offset will suffer. The length of delay used in our experiment is 1 km, or 5 /is. Because of this 
delay, any frequency fluctuation of the LIMO will cause a voltage fluctuation at the output of the mixer. 
We measured the spectrum of this voltage fluctuation with a high dynamic-range spectrum analyzer and 
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transferred the spectral data to a computer. Finally, we converted this information into the phase noise 
spectrum of the LIMO according to the procedures given in [28]. In these experiments, the noise figure 
of the RF amplifier was 7 dB. 

C. Phase Noise as a Function of Offset Frequency and Loop Delay 

Figure 7(a) is the log versus log scale plot of the measured phase noise as a function of the frequency 
offset, /'. Each curve corresponds to a different loop delay time. The corresponding loop delays for curves 
1-5 are listed adjacent to each curve, and the corresponding oscillation powers are 16.33, 16, 15.67, 15.67, 
and 13.33 dBm, respectively. Curve fitting yields the following phase noise relations as a function of 
frequency offset /': -28.7 - 201og(/'), -34.84 - 201og(/'), -38.14 — 201og(/'), -40.61 - 201og(/'), and 
-50.45 - 20 log(/'). Clearly, the phase noise has a 20-dB-per-decade dependence with the frequency 
offset, in excellent agreement with the theoretical prediction of Eq. (31b). 

Figure 7(b) is the measured phase noise at 30 kHz from the center frequency as a function of loop delay 
time, with data points extracted from curves 1-5 of Fig. 7(a) and corrected to account for oscillation power 
differences. Because the loop delay is increased by adding more fiber segments, the open-loop gains of the 
oscillator with longer loops decrease as more segments are connected, causing the corresponding oscilla- 
tion power to decrease. From the results of Fig. 8, discussed below, the phase noise of the LIMO decreases 




Fig. 7. Single side-band phase noise of a LIMO measured at 800 MHz: (a) phase noise 
spectra at different loop delays and their fits to Eq. (31 b) and (b) phase noise at a 30-kHz 
offset from the center frequency as a function of loop delay time. 
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linearly with the oscillation power. To extrapolate the dependence of the phase noise on the loop delay 
only from Fig. 7(a), each data point in Fig. 7(b) is calibrated using the linear dependence of Fig. 8, while 
keeping the oscillation power for all data points at 16.33 mW. Again, the experimental data agree well 
with the theoretical prediction. 

D. Phase Noise as a Function of Oscillation Power 

We have also measured the phase noise spectrum of the LIMO as a function of oscillation power, with 
the results shown in Fig. 8. In that experiment, the loop delay of the LIMO was 0.06 /is, the noise figure 
of the RF amplifier was 7 dB, and the oscillation power was varied by changing the photocurrent, I v h, 
according to Eqs. (3), (4), (18a), (18b), and (18c). With this amplifier and the photocurrent level (1.8 mA 
- 2.7 mA), the thermal noise in the oscillator dominates. Recall that, in Eqs. (25), (26), and (27), the 
phase noise of a LIMO is shown to be inversely proportional to the oscillation power. This is true if the 
gain of the amplifier is kept constant and the photocurrent is low enough to ensure that the thermal noise 
is the dominant noise term. In Fig. 8(a), each curve is the measurement data of the phase noise spectrum 
corresponding to an oscillating power, and the curves in Fig. 8(b) are the fits of the data to Eq. (31b). 
Figure 8(c) is the phase noise of the LIMO at 10 kHz as a function of the oscillation power, extracted 
from the data of Fig. 8(b). The resulting linear dependence agrees well with the theoretical prediction of 
Eq. (31b). 

E. Phase Noise Independence of Oscillation Frequency 

To confirm our prediction that the phase noise of the LIMO is independent of the oscillation frequency, 
we measured the phase noise spectrum as a function of the oscillation frequency, and the result is shown 
in Fig. 9(a). In the experiment, we kept the loop length at 0.28 /zs and varied the oscillation frequency 
by changing the RF filter in the loop. The frequency was fine tuned using an RF line stretcher. It is 
evident from Fig. 9(a) that all phase noise curves at frequencies 100, 300, 700, and 800 MHz overlap with 
one another, indicating a good agreement with the theory. Figure 9(b) is a plot of the phase noise data 
at 10 kHz as a function of the frequency. As predicted, it is a flat line, in contrast with the case when 
a frequency multiplier is used to obtain higher frequencies. This result is significant because it confirms 
that the LIMO can be used to generate high-frequency signals up to 75 GHz with a much lower phase 
noise than that which can be attained with frequency-multiplying techniques. 


V. Summary 

We have introduced a high-frequency, high-stability, high spectral purity, widely tunable electro-optic 
oscillator, which we have termed a LIMO. The high stability and spectral purity of the LIMO result from 
the extremely low energy storage loss realization obtained with a long optical fiber. The optical fiber 
is also virtually free of any frequency-dependent loss, resulting in the same long storage time and high 
spectral purity signals for both low- and high-frequency oscillation. On the other hand, the oscillation 
frequency of the LIMO is limited only by the speed of the modulator, which at the present can be as 
high as 75 GHz. As yet another unique feature, the output of the LIMO may be obtained directly as 
microwave signals or as intensity modulations on an optical carrier for easy interface with optical systems. 

We also analyzed the performance of this oscillator by deriving expressions for the oscillation threshold, 
Eq. (5); oscillation amplitude, Eqs. (18a), (18b), and (18c); and oscillation frequency, Eqs. (19a) and (19b). 
These results agree quite well with experimental data obtained with laboratory versions of the LIMO. 

We also derived the expression for the spectrum of the output of the LIMO and showed that it has a 
Lorenzian line shape, given by Eqs. (25) and (27). The spectral width of the output signal was found to 
be inversely proportional to the square of the loop delay time, given by Eqs. (28a) and (28b). In addition, 
at low optical pumping levels where thermal noise dominates, the spectral line width was found to be 
inversely proportional to the oscillation power of the oscillator, similar to the Schawlow -Townes formula 
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Fig. 9. Single side-band phase noise measurements of the LIMO at different oscillation 
frequencies: (a) phase noise spectra and (b) phase noise at a 10-kHz offset frequency as 
a function of osciiiation frequency, extracted from (a). The loop delay for the 
measurements is 0.28 \xs. 


describing the spectral width of a laser. Since increasing the optical pump power increases the oscillation 
power, in this regime the line width of the LIMO decreases as the optical pump power increases. On the 
other hand, at high pump powers where the pump laser's relative intensity noise dominates, the spectral 
width approaches a minimum value determined by the laser’s RIN noise, as given by Eqs. (25), (27), and 
(34). We measured the phase noise spectrum of the LIMO and verified our theoretical findings. 


It is important to note that the analysis performed here was for the specific case of the LIMO with 
a Mach Zehnder electro-optic modulator. Other modulation schemes, such as with electro-absorptive 
modulators or by direct modulation of semiconductor lasers [29], will also lead to signals with charac- 
teristics similar to those obtained in this work. For these examples, the theoretical approach developed 
above is still applicable after suitable modifications. The major change required in the analysis is the 
replacement of Eq. (1), which describes the transmission characteristics of a Mach-Zehnder modulator, 
with the appropriate equation for the specific modulation scheme. All other equations can then be derived 
in the same way as described in the theory. 
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Because of its unique properties, the LIMO may be used in a number of applications. As a voltage- 
controlled oscillator (VCO) [12], it can perform all the VCO functions in both electronic and photonic 
applications, including generating, tracking, and cleaning RF carriers. The LIMO has the unique property 
of actually amplifying injected signals [12] and, thus, may be used in high-frequency carrier regeneration 
and signal amplification. Other important potential applications of the LIMO include high-speed clock 
recovery [30,31], comb and pulse generation [12], high-gain frequency multiplication, and photonic signal 
upconversion and downconversion [32] in photonic RF systems [33,34]. 
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